DIRECT DARK MATTER EVENT RATES WITH A VELOCITY 
DISTRIBUTION IN THE EDDINGTON APPROACH 
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Exotic dark matter together with the vacuum energy (associated with the cosmological constant) 
seem to dominate the Universe. Thus its direct detection is central to particle physics and cosmology. 
Supersymmetry provides a natural dark matter candidate, the lightest supersymmetric particle 
(LSP). One essential ingredient in obtaining the direct detection rates is the density and the velocity 
distribution of the LSP in our vicinity. In the present paper we study simultaneously density profiles 
and velocity distributions in the context of the Eddington approach. In such an approach, unlike 
the commonly assumed Maxwell-Boltzmann (M-B) distribution, the upper bound of the velocity 
arises naturally from the potential. 

PACS numbers: 95.35. +d, 12.60.Jv 



I. INTRODUCTION 



The combined MAXIMA-1 Q, BOOMERANG 0, DASI and COBE/DMR Cosmic Microwave 
Background (CMB) observations imply that the Universe is flat @, fi = 1.11 ± 0.07 and that 
most of the matter in the Universe is Dark 0. i.e. exotic. Combining the WMAP data with other 
■ experiments, crudely speaking one finds: 

q ; n b = 0.05, ricDM = 0.30, n A = 0.65 

• Since the non exotic component cannot exceed 40% of the CDM pfl, there is room for the exotic 
WIMP's (Weakly Interacting Massive Particles). Supersymmetry naturally provides candidates for 
the dark matter constituents 0-13 • In the most favored scenario of supersymmetry the LSP (Light- 
est Supersymmetric Particle) can be simply described as a Majorana fermion, a linear combination 
of the neutral components of the gauginos and higgsinos 0-0] ■ ^ n m °st calculations the neutralino 
O is assumed to be primarily a gaugino, usually a bino. Even though there exists firm indirect ev- 
idence for a halo of dark matter in galaxies from the observed rotational curves, it is essential to 
directly detect 0-^]| such matter. Until dark matter is actually detected, we shall not be able to 
exclude the possibility that the rotation curves result from a modification of the laws of nature as 
C$ • we currently view them. This makes it imperative that we invest a maximum effort in attempting 
to detect dark matter whenever it is possible. Furthermore such a direct detection will also unravel 
the nature of the constituents of dark matter. 

The possibility of such detection, however, depends on the nature of the dark matter constituents 
(WIMPs). Since the WIMP is expected to be very massive, m x > 30GeU, and extremely non rel- 
ativistic with average kinetic energy T « 50i^eU(m x /100GeU), it can be directly detected 0-0 
mainly via the recoiling of a nucleus (A,Z) in elastic scattering. The event rate for such a process 
can be computed from the following ingredients: 

1. An effective Lagrangian at the elementary par ticle (quark) level obtained in the framework of 
supersymmetry as described , e.g., in Refs (lClll^ . 

2. A well defined procedure for transforming the amplitude obtained using the previous effective 
Lagrangian from the quark to the nucleon level, i.e. a quark model for the nucleon. This step is 
not trivial, since the obtained results depend crucially on the content of the nucleon in quarks 
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other than u and d. This is p articularly true for the scalar couplings, which are proportional 
to the quark masses — Q , as well as the isoscalar axial coupling 0, 0] . 

3. Knowledge of the relevant nuclear matrix elements 0| — E3> obtained with as reliable as 
possible many body nuclear wave functions. Fortunately in the case of the scalar coupling, 
which is viewed as the most important, the situation is a bit simpler, since then one needs 
only the nuclear form factor. 

4. Knowledge of the WIMP density in our vicinity and its velocity distribution. Since the essential 
input here comes from the rotational curves, dark matter candidates other than the LSP 
(neutralino) are also characterized by similar parameters. 

In the past various velocity distributions have been considered. The one most used is the isother- 
mal Maxwell-Boltzmann velocity distribution with < v 2 >= (3/2)uq where vq is the velocity of the 
sun around the galaxy, i.e. 220 km/s. Extensions of this M-B distribution were also considered, in 
particular those that were axially symmetric with enhanced dispersion in the galactocentric direction 
|19ll20| . In such distributions an upper cutoff v esc = 2.84t>o was introduced by hand. 

Non isothermal models have also been considered. Am ong: those one should mention the late infall 
of dark matter into the galaxy, i.e caustic rings 0, 122I I23L EiL |2]| , dark matter orbiting the Sun 
pfij . Sagittarius dark matter |27j . 

The correct approach in our view is to consider the Eddington proposal j2^| , i.e. to obtain both the 
density and the velocity distribution from a mass distribution, which depends both on the velocity 
and the gravitational potential. Our motivation in using Eddington 28] approach to describing the 
density of dark matter is found, of course, in his success in describing the density of stars in globular 
clusters. Since this approach adequately describes the distribution of stars in a globular cluster in 
which the main interaction is gravitational and because of its generality , we see no reason why such 
an approach should not be applicable to dark matter that also interact gravitationally. It should 
be noted that the attempt to use Maxwellian (M-B) distribution to describe the star distribution 
in globular clusters led to results that did not correspond to observations (2^ ■ So it seems that the 
use of M-B distribution to describe dark matter is not very well motivated and a different approach 
is required. 

It seems, therefore, not surprising that this approach has been used by Merritt [2^ and applied 
to dark matter by Ullio and Kamionkowski |30j and more recently by us |3l| . It is the purpose of 
the present paper to obtain a dark matter velocity distribution, which is consistent with assumed 
halo matter distributions and has a natural upper velocity cut off. The results presented here are 
motivated by the dark matter candidate provided by supersymmetry, namely the LSP (neutralino). 
They can easily be extended, however, to be applied to other heavy WIMP candidates. 

II. THE DARK MATTER DISTRIBUTION IN THE CONTEXT OF THE EDDINGTON 

APPROACH 

As we have seen in the introduction the matter distribution can be given as follows 

dM = 2n f(Jf>(r),v r ,Vt) dx dy dz vt dvt dv r (1) 

where the function / the distribution function, which depends on r through the potential $(r) and 
the tangential and radial velocities vt and v r . In general the distribution function is not symmetric. 
In the above expression we assumed that it is only axially symmetric, with the two tangential 
components being equal. Thus the density of matter p satisfies the equation: 

dp = 27r /($(r),iv,u t ) v t dv t dv r (2) 

It is more convenient instead of the velocities to use the total energy E and the angular momentum 
J via the equations 



J = v t r , 2E = vl + — + 2 $(r) 



(3) 



The use of these variables, which are constants of motion, is very useful, when one wants to study 
steady states. In doing this Eddington used the result of Jeans js3 that the density must be a 
function of first integrals of the equations of motion which follows from Liouville's theorem. The 
advantage of this approach is that the density can be 'inverted' and the velocity distribution can be 
found. Following this approach we find 
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f(E,J) J 



y/2(E-Q(r))- J 2 /r 2 



dJ dE 



(4) 



The limits of integration for E are from $ to and for J from to [2r 2 (E — $(r))] 1//2 . 
Furthermore if the distribution function is known one can obtain the velocity distribution at some 
point, e.g. in our vicinity, by 

/($(r),U r ,U t )|r=r s 

The problem which is more interesting is: Can one obtain the distribution function given the density 
(and hence the potential via Poisson's equation)? The answer is affirmative via Eddington's treat- 
ment of the distribution function and quite easy, if the distribution does not explicitly depend on 
J, but is only a function of E. In the present work we will be concerned with spherically symmetric 
velocity distributions and we will leave out the more realistic axially symmetric case |33j.|34j. Such 
axially symmetric velocity distributions have, however, been found to have interestin g co nsequences 
on the direct dark matter detection rates, especially in directional experiments |35|. |36|. 
If the angular momentum dependence is ignored, by integrating Eq. one finds 



Air J f{E)^2(E-<£(r))dE 



(5) 



with the range of E as above. In this case one can obtain the density as a function of the potential. 
Conversely if the density is given as a function of the potential one can proceed to find the distribution 
function according to the Eddington approach. The distribution then is a function of the total 
energy E = v 2 /2 + <E>(r) and satisfies the Boltzmann Equation with the collision term zero, i.e. 



(V- Vr - V*-Vv)/ = 

In this case the distribution can be expressed as follows: 

d<S> 

IE V* 3 

The above equation can be rewritten as: 
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(6) 



(7) 



(8) 



Thus we can obtain the distribution, if the density p is given as a function of the potential $. We 
find it convenient to rewrite the last equation in terms of dimensionless variables by introducing: 



p_ 

Po' 



$0 



(9) 



The constants po and $0, which set the scale of these two quantities, are related through Poisson's 
equation. Then the Eddington distribution function takes the form: 
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(10) 



where e is given by 



i.e. e the negative of the total energy e (e > 0. In the presence of asymmetry one finds 

1 V + a s u f 2 ) (12) 



* 2|#o| 

where a s is the asymmetry parameter, vt is the tangential velocity. 

For numerical integrations it is more convenient to rewrite the last integral as follows: 



f( \ Po 
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(13) 



The first term is singular as e — ■> 0. This, however, causes no problem, since, as we have seen in 
our earlier work, the integrals over the velocity distribution are relevant in dark matter calculations 
and these remain finite as the velocity approaches the maximum velocity. Anyway in the example 
considered in the present work 77 (0) = 77 (0) = 0. 

Once this function is known one obtains the velocity distribution of matter in our neighborhood 
(r = r s = 0.8a , a = the galactic radius)with respect to the center of the galaxy via the relation 

fv(v) =/(e)| r=r , =0 .8o J (14) 

which must be normalized. The characteristic feature of this approach is that the velocity distribu- 
tion vanishes outside a given region specified by a cut off velocity v m , by setting e| r=rs= o.8a = 



III. A SIMPLE REALISTIC DENSITY PROFILE 

We will consider three types of matter density: 

• A spherical ordinary matter density (bulge density) 

• Ordinary matter density in the form of a disc 

• Dark matter density. 



A. Spherical ordinary matter density 

To obtain analytical expressions we will simulate this density as follows: 

M*)=P0b (1 + J )5/2 , x= T - (15) 

with pob a constant and a the galactic radius. So if this were the whole story, p b — 2p s with 
p s the mass density in our vicinity. A distribution found by Plummer (3lJ and Von Zeipel [38|. 
p{x) oc (1 + a; 2 ) -5 / 2 was obtained by considering a gas in spherical container with a specific heat 
ratio of 7 = 1.2 . We note that a value of 7 = 1, which leads to isothermal distribution, was excluded 
by observation. This normal matter density profile is shown in Fig. ^ With this density one finds 
the potential: 

* h (x) = -4TrG N a 2 Pob - ^— (16) 
ova: + 1 



and the rotational velocity: 



V " {X) = X — = AlrGNa 3 (^2 + 1)3/2 



The above potential and rotational velocity are shown in Fig. 




FIG. 2: The Potential and the rotational velocity in units of 47rGjva 2 pof> due to the spherical ordinary matter 
distribution discussed in the text. 



B. Ordinary matter density distributed on a disc 

This is a more complicated problem. We will adopt the simplification that the matter distribution 
is a S function along the axis of the disk. The radial density has been successfully modelled in the 
form of exponential profiles [39j . Nevertheless, to minimize the number of parameters employed, we 
will assume further that it has the same radial dependence as discussed above. In other words 

y/2 r 

Pd(x,z) = podS(z)— — —^j- , x = ~ (18) 
(1 + x z ) t >* z a 

where r is now the radial distance from the axis of symmetry (both x and z are measured in units 
of a). The potential now takes the form: 

1 f°° 

§d{x,z) — AttGno 2 pod— / dk cos (kz)G(x,k) (19) 

7T Jo 

where the Green's function in momentum space is given in terms of the modified Bcsscl functions: 

G(x,k) = G<(x,k) + G > (x,k) (20) 

with 

G > = -I (kx) K (ky) Pb (y)dy (21) 

J X 

G < = -K Q (kx) f I (ky)p b (y)dy (22) 
Jo 



The radial rotational velocity is given by 



1 



v d {x, z) — A-kGno- Pod—x / cos(kz)g(x, k)dk 



(23) 



with 



g(x, k) = g < (x, k) + g> (x, k) 



(24) 



9< = Io(kx)K (kx)p d (x) + kKi(kx) / I (ky)p d (y) dy (25) 

Jo 



g< = I (kx)K (kx)p d - kh(kx) \ K (ky)p d (y)dy 



(26) 



The obtained potential and the square of the rotational velocity on the plane of the galactic plane are 
shown in Fig. |3J In the same Figure we also plot the same quantity obtained with the exponential 
profile: 



p ed (x,z) = p 0ed 5(z)e 



-3.5x 



r 

X = — 

a 



(27) 



We see that the two densities give essentially the same rotational velocities with the possible excep- 
tion at very small distances. 




FIG. 3: On top we show the Potential in units of 4-kGno, pod on the plane of the galactic disc resulting 
from the ordinary matter distribution on the disc discussed in the text. At the bottom we show on the left 
the square of the rotational velocity resulting from the above density and on the right the same quantity 
resulting from an exponential density profile. 



C. Dark matter density 



There are many halo density profiles, which have been employed. In the present case we will 
consider only spherical distributions, since it is not easy to extend the Eddington approach to deal 
with the most general case. Among the most commonly used are: 



The simple density profile 



1 + ar a 



with a the radius of the Galaxy. This profile has the advantage that the rotational velocity 
remains constant with the distance from the center of the galaxy becomes very large. 



Another simple profile is: 



p{x)= ,/° x=- (29) 
' K ' x(l + x) 2 a v ; 



suggested by N-body simulations [2(j] • This profile provides a better description of the expected 
density near the center of the galaxy. It does not, however, predict the constancy of the 
rotational velocities at large distances. 

In the present work we will consider the density profile of Eq. (|28|l 

Unfortunately with this density the potential diverges at infinity. On the other hand the solution to 
Poisson's equation is finite at the origin and it can be chosen to vanish there, i.e. it takes the form: 

— = ^ + -lnl+^-1 30 

q?o x 2 

One may choose a radius x — c outside of which the density can be chosen to go faster to zero. One 
convenient choice is: 

with the requirement that at x — c the density is continuous with a continuous derivative, we thus 
find: 



P>(x) = Po 



2(c 2 + l) (c 2 + l) 2 



(32) 



(x 2 + iy {x 2 + 1) 

The obtained density is shown in Fig. 0] The potential in the region x -< c is the same as before: 

= /<(*)+<* (33) 
The solution in the outer region takes the form: 

%M = 7< (c ) + 7> {x) _ 7> (c ) _ C3 (34) 

where the constant C3 can be chosen to make the potential vanish at some point and 

tan _1 (a;) 1, , a ,„_ s 

J< = ^ + - log (x 2 + 1) 35 

x 2 



l/(c 2 + l) 2 (c 2 -7)tan- 1 (x)(c 2 + l) (-c 4 + 6c 2 + 15) tan^c) - c (c 2 + 15) ' 



(36) 




FIG. 5: The potential £ = 5£^, 'i.e the absolute value of the potential in units A-kGno, 2 po, obtained with 
the density function shown in Pig. 0]and chosen to vanish at infinity. The potential drops to zero very slowly 
at large distances. 



by choosing C3 as 

C3 = i(21og( C 2 + l)+7) 

the potential can be made to vanish at infinity. It is shown in Fig. With the above ingredients 
it is not very difficult to obtain the function needed in the Eddington approach, namely 77(C)- The 
results are shown in Fig. In all cases Poisson's equation yields a relation between $0 and po, 
namely 

$0 = — 47rGjva 2 po- 



The obtained rotational velocity curve, due to dark mater alone, is given in Fig. m units of 
\J AttGno^pq From the observed rotational velocity one can fix the constant po. 
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FIG. 6: The density r\ as a function of the potential £ 




FIG. 7: The rotational velocity due to dark matter as a function of the distance in units of v/ AtvGnci 2 po- 
Shown on the left is the one obtained with the density profile of Eq. 1281 adopted in this work, while on the 
right the NFW profile see Eq. was employed. 



D. Combining ordinary matter and dark matter 

We will now combine the three kinds of distribution considered above. We will assume that 
Pod = Pob = Ps an d p s — Po- This means that in our vicinity the density of dark matter is equal 
to that of ordinary matter. Furthermore the ordinary matter density in our vicinity is equally split 
between the spherical and disc geometries. The obtained results are shown in Fig. 18181 We see that 




FIG. 8: The Potential due to all three matter distributions considered in the text. 



in order to fit the rotational velocity of the sun, 220 km/ s we need a density po = 1-0 x 10~ 21 kg/m 3 . 
This is in agreement with the value of 0.3 GeV/cm 3 = 0.5 x 10 -21 kg/m 3 used in calculations of 
the rates for direct dark matter searches (in our model the dark component is half of the total). 




FIG. 9: The rotational velocity resulting from the potential of Fig. |H| 



IV. THE VELOCITY DISTRIBUTION IN OUR VICINITY 



It is clear that the velocity distribution and, in particular the maximum allowed velocity, is related 
to the escape velocity via the density p . For the moment we will ignore the asymmetric term and 
set a s = 0. 



A. Dark Matter only 

The relation between the density and the potential has already been shown in Fig. H3 In this case 
we obtain the velocity distribution shown in Fig. 1101 with the usual normalization imposed 



y 2 fv(y)dy = 1. 



Instead of the velocity we have used the dimcnsionlcss quantity y: 



y/4irG N a 2 p 

For comparison we present the velocity distributions obtained with the profile of Eq. I|29[) in Fig 




y — * 



FIG. 10: The velocity distributions f v (y) and y 2 fv(y) in units of the parameter yj 4-kGno? pa, i.e. y = 
vj \l AkGncl 2 po obtained with the density profile of Eq. I28H 



ITT1 and with the standard M-B distribution in Fig. [T^l The maximum velocity allowed by our 
distribution is only y m = 2.8. Assuming a = 3.1 x 10 20 m and p\ x =i — Ps/2, i-e. half of the 
local density to be due to dark matter, we find p = p s which yields v m = y m ^AttGn^Po = 
2.8 x 270km/ s = 7.5 x 10 5 m/s. This is a bit higher than the escape velocity, v esc = 6.2 x 10 5 m/s, 
assumed in theories employing the M-B distribution. For comparison we present the same quantities 
for the M-B distribution in Figs ll2|j"2l Since the M-B distribution does not go to zero at finite values 



FIG. 11: The velocity distributions f v (y) and y 2 f v {y) in units of the parameter -J A-kGncl 2 po, i.e. y = 
vj \J 'A-kGncl 2 po, obtained with the density profile of Eq. H29H 




y — * 

FIG. 12: The same as in Fig. ESI in the case of the M-B distribution. Now, however, the parameter y is the 
velocity in units of the sun's rotational velocity, i.e. y — v/vo with vo = 2.2 x 10 m/s. Note that, in these 
units, the escape velocity is at 2.84 



of the velocity, the maximum allowed velocity is set by hand equal to the escape velocity. 

We should emphasize that in this approach we encounter two characteristic velocities. One is the 

rotational velocity and the other the maximum allowed velocity v m . The first depends on the square 

root of the derivative of the potential, while the second scales with the square root of the potential 

itself. 



B. Both ordinary matter and dark matter 



To improve the situation one will attempt to include gravity. Since it is very hard to incorporate 
the disc geometry into the Eddington approach, we will attempt to mimic the gravitational effects 
of ordinary matter with a spherical distribution like the one discussed above, but twice a large, so 
that in our position the usual and cold dark matter contribution are about equal. The thus obtained 
function 77 = r](£) is shown in Fig. 1131 

In this case we find y m = 3.2. In other words the effect of ordinary matter is small, since the 
dark matter potential in our vicinity is about 12 times stronger than the potential due to ordinary 
matter. Again the condition p\ x =i — ps implies that overall constant po m the density distribution 
is p s , which gives v m = 3.1 x 2.7 x 10 5 m/s = 8.0 x 10 5 m/s. With this modification the obtained 
results for the velocity distribution are shown in Figs [21 We see that the inclusion of gravity has 
very little effect on the velocity distribution. So in what follows we will consider only the dark 
matter component. 
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FIG. 13: The density 77 as a function of the potential £, when both ordinary and dark matter are included. 
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FIG. 14: The same as in Fig. I1UI when both ordinary and dark matter are included. 



V. DIRECT DARK MATTER RATES 

In the present work we find it convenient to write the event rate in the form 1 : 



R=K 



C co h(A^ r {A))a p x o + C sp in(A, fJtr(A))(T^o Csp 



(37) 



In this expression afj x0 is the LSP-nucleon scalar cross section, while cr^™ is the proton cross 
section associated with the spin. The quantity Cspmis given by [l5| : 



1 



3(1 + 4) 

J A 



— S(u) 



(38) 



s(u) » 5(0) = [An Q (o)) 2 + 2^r» (o)fi 1 (o) + r»!(o)) 2 ; 



n 



(39) 



The static spin matrix elements are obtained in the context of a given nuclear model (see, e.g., 
previous work 01 — and references therein). Even, though most of what we are going to 

say applies in the case of the spin induced rate we are not going to further elaborate here. 



see our previous recent work ll r ll4(l for a more detailed discussion and additional references 



In Eq. (S3 K is given by 
K 

and 



K7^— V(^>-160 10 ( P b) y 



100 GeV m E 



c co/l (A,/i r (A)) = 



100 GcV 



/x r (A) 
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Mr(p) 



v4 t coh (A) 



(40) 



(41) 



lOOGeV 



Hr{p) 



A 



(42) 



The parameters c co h(A, fi r (A)), c sp i n (A, fi r (A)), which give the relative merit for the coherent and 
the spin contributions in the case of a nuclear target compared to those of the proton, have already 
been tabulated jljj f° r energy cutoff Q m in = and 10 keV. 

Via Eq. I|37l) we can extract the nucleon cross section from the data. The most interesting quantity 
is t coh {A). It is defined as: 



tr 



dt r 



du 



-du 



(43) 



u is the energy transfer to the nucleus (in dimensionless units, see below) 

u m in & detector threshold 



Umax maximum WIMP velocity 



dt co h _ 2 V par . 2| eV \|2 



du V 3 Do 

for the coherent mode and 



T(u) , T(u) = a 2 \F(u)\ 2 V(aVu) (44) 



dt^n = /2Vr T(tl) T{u) = ^p^^ (46) 

du V o wo 

where v par — v m in the present approach and v par = vq in the case of the M-B distribution and The 
nucleon cross sections, which carry the dependence on the particle model parameters, are the most 
important ones, but they are not of interest in our present calculation. One such parameter is, of 
course, the WIMP mass. In the above expressions F(u) is the form factor, entering the coherent 
scattering and F±±(u) is the spin response function entering via the axial current. The function \1/ 
depends on the WIMP distribution velocity employed and is a function of the energy Q transferred 
to the nucleus 

U =7T > Q° = = 4- 1 x 104 ^ 4/3 KeV ( 47 ) 

where A is the nuclear mass number and the dimensionless parameter a is given by: 



a = [V2^L r bv par ] 



(48) 



where \x T is the reduced mass of the WIMP-nucleus system and b is the (harmonic oscillator) size 
parameter. 

The function, which is basic to us, 'S, is given by 

*(aV^)= dy [ sin Odd f d<f> y /„(y,M) (49) 

Jay/u Jo Jo 

with 

f v (y,6,(f>) = f v (y/y 2 + 2yy sun cos6 + y^ un ) (50) 

where 9 is the polar angle as measured from the direction of the sun's motion and y sun is the sun's 
velocity in units of v m . Note that since the argument of the function /„ is constrained to be less 
than y m . If y m is small the allowed region in the (y, £) space is very restricted. The function ^(a^/u) 
is plotted in Fig. ^| The function \& for both distributions is shown in Fig. ^| We should stress 




FIG. 15: On the top we show the function ^(a^/u) which enters the differential (with respect to the energy 
transfer u) event rate in dark matter searches. It has been obtained in the context of the Eddington approach 
as discussed in the text. For comparison we present at the bottom the same function (continuous curve) 
together with that obtained using the M-B distribution (dotted curve). 

that the relative differential rate Jj^j can be obtained by combining the above results of 9(ay/u 
with the nuclear form factor for the target of interest. The dependence on the WIMP mass comes 
through the parameter a. 

VI. APPLICATION IN THE CASE OF THE TARGET 127 I 

In this section we are going to apply the above formalism in the case of a popular target, 127 I, 
which is an odd mass target and can detect both the coherent and the spin modes of the WIMP- 
nuclear interaction. We will include in our only the coherent mode but we do not expect any real 
differences as far as the quantity i ^ is concerned. The nuclear form factor employed was obtained 
in the shell model description of the target and is shown in Fig. The quantity dt *° h obtained 
both for our distribution as well as for the familiar M-B distribution is shown in Figs I17I2UI for 
various WIMP masses m x . It is clear that the differential rate is a fast decreasing function of the 
energy transfer. This is particularly true for low WIMP masses. 






FIG. 17: On the left we show the quantity 



keV 



for ttu 



Q n 



64 



10 GeV in the case of the distribution obtained 



in this work. On the right we show the same quantity in the case of the M-B distribution. The expression 
for tspin is similar. Note that for such a small WIMP mass the differential rate drops very fast as a function 
of the energy transfer. 
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FIG. 18: The same as in Fig. [T7|fbr m x = 30 GeV (solid line) and m x = 50 GeV (dashed line). 



Integrating the differential rate from zero to u max = Umax/a 2 , wltn Umax = Vm {y m ax=2.$,i) for 
the present (M-B) distributions respectively, we find the total rate | as a function of the WIMP 
mass. The results are shown in Fig. [21 As we have already mentioned the nuclcon cross sections 
also depend on the LSP mass. So the absolute event rates, which include these cross sections, are 
expected to drop even faster as a function of the LSP mass. In practice, however, the detectors 
have a low energy threshold. So only the event rates above an energy transfer Q t h can be detected. 
Thus we present the relative total rates t co h as a function of Q t h in figs 1221241 From these plots 
we see that it is crucial for experiments to lower the threshold energy as much as possible. This is 
particularly true for small WIMP masses. 
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FIG. 19: The same as in Fig. E|for m x = 75 GeV (thick solid line), m x = 100 GeV (thin solid line), 
m x = 125 GeV (short dashed line) and m x = 150 GeV (long dashed line). 



0.2 0.4 0.6 0. 




u = Qo = 64 



keV 



FIG. 20: The same as in Fig. EH for m x = 200 GeV (thick solid line), m x = 250 GeV (thin solid line), 
m x = 350 GeV (short dashed line) and m x = 500 GeV (long dashed line) . 




FIG. 21: On the left we show the quantity t co h in the case of the present distribution as a function of the 
WIMP mass and threshold energy {Qth = <=>■ thick solid curve, Qth = 5 keV fine solid curve and 
Qth = 10 keV dashed curve). On the right we show the same quantity in the case of the M-B. It is 
clear that the rates decrease as the threshold energy increases. This is is especially true for low LSP mass 




FIG. 22: The relative rates t co h as a function of threshold energy for WIMP masses m x = 10 (thick solid 
line), m x = 30 GeV (fine solid line) and m x = 50 GeV (dashed line). The results on the left correspond to 
the present distribution, while those on the right to the M-B distribution. 
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FIG. 24: The same as in Fig. [331 for m x = 200 GeV (thick solid line), m x = 250 GeV (thin solid line), 
m x = 350 GeV (short dashed line) and m x = 500 GeV (long dashed line). 



VII. MODULATION 



In the above discussion we did not take into account the motion of the Earth. The expected event 
rates, however, are very small due to the smallness of the nucleon cross sections not discussed in this 
work. So the experiments must fight against formidable backgrounds. Fortunately there are some 
signatures of the WIMP-nuclear interaction, which must be exploited. One such comes from the 
fact that the event rates depend on the relative velocity between the WIMP and the target. The 
most important velocity dependent contribution comes from the rotation of the Earth around the 
sun with a velocitywi = 0.27 vq. It turns out that only the component of the Earth's velocity along 
the sun's direction of motion, {v\) z = 0.135 cosa, is relevant (a is the phase of the Earth, a = 
around June 3nd). We can thus apply the above formalism by 

Vsun 2/«m(l + 0.068 cos a) 

Thus Ea s PSI and ETUI become 

T(u) =S> a 2 \F{u)\ 2 [V(ay/u) + H(aVu)cosa] (51) 

for the coherent mode and similarly for the spin We will examine the modulation effect in the case 
of 127 I. The behavior of the function H(a^/u) is exhibited in Fig. [23 We notice the sign change of 
H as the energy transfer changes. Integrating the differential rates we obtain the total rates, which 
now now take the form: 

t =>• t(l + ft, cos a) (52) 

of course, on the parameter a. Quite generally for a light nuclear system and relatively heavy WIMP, 
the parameter a is large, i.e. the lower part of the Fig The exact behavior of the modulation of event 
rates depends, of course, on the parameter a and the nuclear form factor. Quite generally for light 
nuclear targets, i.e. large a, the lower range of Fig. [22] does not enter and the modulation h is not 
suppressed. On the other hand, for intermediate and heavy nuclei, the lower part tends to cancel 
the effect of the upper part of Fig. [23 This leads to a reduction of the modulation of the total rates 




FIG. 25: On the left we show the function H(ay/u) governing the differential modulation amplitude as a 
function the energy transfer u. On the right we show the unimportant higher order effects ~ cos 2a. We 
note the sign change as the energy transfer increases. 



and even in a change of sign of h, i.e. minimum in June and maximum in December. The behavior 
of the quantity h, in the case of the present velocity distribution foe 127 1, is exhibited in Fig. We 
see that the modulation predicted by the present distribution is much smaller than that obtained 
with a M-B distribution (see Fig. 127( 1 . This is true regardless of the energy cut off. Similar results 
have been obtained in a modified M.B. distribution [4l| obtained in a model that couples gravity 
to a scalar field. This way the characteristic velocity is not vo, the sun's rotational velocity but 
-n Vd y v2 ^> vq. We should emphasize that h is the ratio of the modulated rate to the average 




GeV 



FIG. 26: The modulation function h is shown as a function of the WIMP mass for various threshold energies 
obtained with the present velocity distribution. In this figure: Qth = thick solid curve, Qth ~ 5 keV 
•<=>■ fine solid curve and Qth = 10 keV <=>■ dashed curve. 



rate. In other words the increase of h as the threshold energy increases, results from the fact that 
the average rate (denominator) decreases much faster than the modulated rate (numerator) . So the 
increase in h, which is a good signal against background, is not cheap. It comes at the expense of 
the number of counts. 




(GeV) 

FIG. 27: The modulation amplitude ft as a function of the WIMP mass obtained with a M-B distribution 
in the case of 127 I for Qmin — on the left and Qmin = 10 keV on the right. For the definitions see text. 



VIII. CONCLUSIONS 



In the present paper we discussed the dependence of the direct dark matter event rates on the 
density and velocity distribution of WIMP's. This was done in a self consistent way by applying 
Eddington's approach. Even though at present only spherically symmetric matter distributions can 
be treated this way, one can draw the following conclusions: 

• Simple dark matter distributions can adequately describe the rotational curves of both spheri- 
cal and disc galaxies. With the assumed form of the density its scale was determined by fitting 
the rotational velocity of the sun. The value of dark matter density extracted is in agreement 
with that hitherto employed in dark matter calculations, ~ 0.3 GeV/cm 3 . 

• From the assumed density profile we obtained the distribution function. Evaluation of the 
distribution function in our vicinity yields the velocity distribution. This velocity distribution 
automatically vanishes for velocities larger than a velocity v m . v m depends on the square root 
of the potential in our vicinity. It is quite different from the rotational velocity, which depends 
on the square root of the derivative of the potential. 

• The velocity distribution obtained is not very different from the standard M-B in the range of 
velocities of interest to dark matter searches. 

• We focused our attention on the factor t entering dark matter rates, which is independent of 
the nucleon cross section, i.e. it holds for all heavy WIMPs. The values obtained in our model 
are larger than those obtained using the M-B distribution. The relative magnitude depends 
on the WIMP mass. For small WIMP masses the obtained results are 3-4 times larger than 
those of the M-B case. 

• We also computed the modulation amplitudes H(u) and h. We find that both are more 
suppressed than those obtained with the M-B distribution. It seems that the precise value 
of the modulation sensitively depend on the assumed velocity distribution (the form of the 
function H{u) is independent of the assumed nuclear form factor). It appears, though, that 
the modulation obtained with the present velocity distribution is perhaps too small to be of 
practical interest experimentally. 

We did not consider in this work asymmetric velocity distributions. Another signature not discussed 
in this work is the asymmetry, with respect to the sun's direction of motion, expected in directional 
experiments |4^ . In such experiments, which detect not only the energy of the recoiling nucleus, but 
its direction of motion as well, the modulation is expected to be direction dependent. It is expected 
to be quite large in some directions. Such effects in the context of the present approach are currently 
under study. 
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